Neutrino signatures of the supernova - gamma ray burst relationship 



Soebur Razzaque,^ Peter Meszaros^ and Eli Waxman^ 
^Dpt Astronomy & Astrophysics, Dpt Physics, Pennsylvania State Univ., University Park, PA 16802, USA 
^Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel 

(Dated: February 2, 2008) 

We calculate the TeV-PeV neutrino fluxes of gamma-ray bursts associated with supernovae, based 
on the observed association between GRB 030329 and supernova SN 2003dh. The neutrino spectral 
flux distributions can test for possible delays between the supernova and the gamma-ray burst events 
down to much shorter timescales than what can be resolved with photons. As an illustrative example, 
we calculate the probability of neutrino induced muon and electron cascade events in a km scale 
under-ice detector at the South Pole, from the GRB 030329. Our calculations demonstrate that km 
scale neutrino telescopes are expected to detect signals that will allow to constrain supernova-GRB 
models. 
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I. INTRODUCTION 

Gamma-ray bursts (GRB) have recently been shown to 
be associated with supernova (SN) events. This is based, 
most notably, on the burst GRB 030329, in whose optical 
afterglow the supernova (labelled SN 2003dh)_photomet- 
ric and spectral signatures show up 9.6 days Q after the 
7-ray trigger. This late brightening of the supernova is 
caused by the large initial opacity to optical photons of 
the envelope of the exploding supernova. The collapse of 
the stellar core is thought to lead, for progenitor stellar 
masses in excess of 30 Mq, to a blackhole responsible for 
the GRB (e.g. for a review), while for less massive 
progenitors it leads to a rotating neutron star, typically 
detectable as a pulsar. A question of great interest is 
whether the supernova explosion of a progenitor in ex- 
cess of ~ 3OM0 can lead to a black hole promptly (on a 
timescale comparable to the core collapse free-fall time) , 
or else after some delay, during which an initial pulsar 
accretes additional gas which falls back onto it from the 
envelope Q • In the extreme version of this scenario the 
collapse to a black hole is delayed by weeks to months 
Q. Optical observations cannot test the time off-set 
which characterizes this scenario to an accuracy better 
than a few days P, 0, 01 j due to observational difficul- 
ties and uncertainties in modeling the diffusive photon 
transport in the dynamically evolving envelope. Neutri- 
nos, however, are not subject to the very long diffusion 
times of photons, due to their much lower opacities, and 
should provide a much tighter time-tracking of the col- 
lapse. While thermal (MeV) neutrinos produced in the 
collapse are presently undetectable except from the very 
nearest galaxies, ultra-high energy (> TeV) neutrinos 
from GRBs are expected to be detectable from cosmolog- 
ical distances with kilometer scale Cherenkov detectors 
[L> A -2j 2M or giant air-shower detector ^fj ■ Here we 
calculate the > TeV neutrino signatures associated with 
the collapse of massive stellar progenitors of GRB, both 
in the case where a supernova envelope is ejected simul- 
taneously with, or sometime before, the GRB event. We 
calculate the resulting muon and electron cascade events 



in kilometer scale sub-ice detctors, and discuss the char- 
acteristic signatures associated with various time off-sets 
between the GRB and the SN. 



II. OBSERVATIONS AND MODEL 



The GRB 030329 is located at the coordinates a — 
10''44'"50.03", S = +21°31'18.1" (J2000) Q. The dis- 
tance of the burst at redshift z « 0.17 is D 10^'^-^D27.3 
cm for an = 0.3, JIa = 0.7 and Ho = 75 km 
s~^ Mpc~^ cosmology. The observed GRB duration is 
At - 30Ai3o s with 10"^ erg/cm^ fluence in the 30- 
400 keV range. The peak flux is 7 x 10~^ erg cm~^ 
s~^ with 1.2 s duration. The isotropic 7-ray luminosity 
L!^° ^ 10^^L-y_5i erg/s is then approximately the same as 
the typical value for long bursts. 

We take the progenitor star to have an original mass 
of ~ 30 — 4OM0, which subsequently lost its H envelope 
leaving behind a He core of ^ 14M0 in its presupernova 
phase. The He core has initially a smaller radius of ~ 
1.2i?0, for a solar metallicity star, but grows larger later 
in its He burning stage. The density profile of the star can 
be roughly estimated as p = p*(i?,/r — 1)", where p* = 2 
g/cm^ and i?* ~ 10^^ cm in the presupernova phase [T^ . 
In case of 71 w 3, for radiative envelope, p « 10"'^ g/cm^ 
for r - 0.9i?*. 

The supernova SN 2003dh is similar to SN 1998bw, 
with a similar explosion energy estimated as E^n = 
(2 — 5) X 10^^ erg, and an ejected supernova remnant 
(SNR) shell velocity of « 0.1c @. We consider here two 
scenarios: first, the SN takes place 0.1-8 days prior to 
the (elcctromagnctically detected) GRB event; and sec- 
ond, both the SN and the GRB occur simultaneously. 
The neutrino fluxes from these two scenarios differ sig- 
nificantly, as discussed in the following section. 
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III. NEUTRINO FLUX COMPONENTS 

High energy neutrinos are created from photomeson 
{pj) and proton-proton (j>p) interactions by the shock 
accelerated protons in the GRB shocks. In the 7- 
ray prompt phase of the GRB, neutrinos are produced 
through p7 interactions (13j . A precursor component 
arises, due to pp and p^ interactions, when the GRB 
jet is burrowing it's way through the star If an 

SNR sheU is ejected days before the GRB event, an extra 
component is added to the prompt phase neutrino flux 
becasue of pp and pj interactions in the shell mill. 
Afterglow reverse shocks also generate high energy neu- 
trinos through p7 interactions 18, 19]. We consider each 
of these flux components, in the framework of astrophys- 
ical models described in Sec. Illl below. 



Burst 



The Thomson depth of the GRB internal shocks at a 
radius rgh is rxh ~ crThU^°6t/{4:Trr^^J^mpc'^), where St ~ 
lO^^St^^ s is the variability time and F w SOOFsoo is the 
bulk Lorentz factor. Solving for rsh = T-thin at tth — 1, 
we get the shock radius at which the GRB jet becomes 
optically thin as 
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The optical depth for p7 — > A at threshold with observed 
gburst peak synchrotron photons is 



(2) 



Correspondingly, the proton to pion conversion efficiency 
in the internal shocks is f^™'^^ ~ If^™^*". The neutrino 
break energy for p^ interactions with peak synchrotron 
photons is 
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assuming 5% of the proton energy goes to a neutrino in 
each interaction at the A-resonance (E'^e!^"''^' = 0.2 GeV^ 
in the comoving frame). Neutrinos below this energy in- 
teract with synchrotron photons above e!^"''''* w 0.5 MeV 



which follow a typical power law distribution dN^/dE. 
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At ultra- high energies, neutrino production is sup- 
pressed due to synchrotron losses of secondary pions and 
muons produced from p^ interactions. The maximum 
pion synchrotron break energy, from equipartition mag- 
netic field with equipartition fraction — 0.01^b,-2, 



Assuming the luminosity of internal shock-accelerated 
protons in the GRB is Lp ~ kL^ 51, the power-law dis- 
tribution of protons is 
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Here ^ 1 is the fraction of protons accelerated and 
K > 1. The corresponding neutrino spectrum at earth 
= d'^N^/dE^dt = {4:T:D^)-^d^Np/ dEpdt] is 
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where E^''^^^ 



9.4 X 10^ GeV is the maximum syn- 
chrotron energy for neutrinos from Eq. Q). 

B. Afterglow 

The GRB afterglow arises as the jet fireball ejecta 
runs into the ambient medium [the inter-stellar medium 
(ISM) or a pre-ejected stellar wind], driving a blast wave 
ahead into it and a reverse shock back into the GRB 
jet ejecta. This (external) reverse shock takes place well 
beyond the internal shocks, at a radius Ve ^ AV^c/^t ^ 
lO^^'^rigpAtao cm Isj. Here « 250r25o is the bulk 
Lorentz factor of the ejecta after the partial energy loss 
incurred in the internal shocks. Neutrinos are produced 
in the external reverse shock due to P7 interactions of 
internal shock accelerated protons predominantly with 
synchrotron soft x-ray photons produced by the reverse 
shock. The characteristic photon energy in the observer 
frame for the ISM case is [l§l 
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where 0.01^b^_2 is the magnetic field equipartition frac- 
tion from typical afterglow fits and Ino cm~'^ is the 
inter-stellar medium density. The corresponding neu- 
trino break energy, for p^ interactions, is 
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The efficiency of pion conversion from p^ interactions 
in GRB afterglows is /|'°" ~ 0.1/^'°^; smaller than in 
internal shocks 'iSj. The afterglow neutrino spectrum 
[Afterglow (ISM) case] is then 
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Note that ultra-high energy neutrinos from afterglows are 
not suppressed due to pion synchrotron losses, because 
of the smaller magnetic fields at the larger radius of the 
external shock (in constrast to internal shocks). 

The estimate in Eq. is derived under the assump- 
tion that the jet expands into a medium with a density 
typical to that of the interstellar medium, rt ~ 1 cm~'^. 
However, in the case of a massive star progenitor the 
jet may be expanding into a wind, emitted by the pro- 
genitor prior to its collapse. In this case, the density of 
the surrounding medium, at the external shock radius, 
may be much higher than that typical to the interstellar 
medium. For a wind with mass loss rate of 10~^Mq yr~^ 
and velocity of = 10^ km/s, the wind density at 
the typical external shock radius would be ~ 10^ cm~'^. 
The higher density implies a lower Lorenz factor of the 
expanding plasma during the reverse shocks stage, and 
hence a larger fraction of proton energy lost to pion pro- 
duction. Protons of energy E„ ^, 10^* eV lose all their 
energy to pion production [l^, [13 > and the resulting neu- 
trino flux [Afterglow (wind) case] is approximately given 
by M 
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The neutrino flux is expected to be strongly suppressed 
at energies E^ > 10^^ eV, since protons are not expected 
to be accelerated to energies Ep :s> 10^° eV. 



C. Supranova 

If the SN explosion resulting in a SNR shell ejection 
takes place hours to few days prior to the GRB ("supra- 
nova" scenario), then one expects an extra neutrino com- 
ponent, due to GRB-accelerated proton interactions tak- 
ing place in the pre-ejected SNR shell, which is in ad- 
dition to the burst and afterglow components. Optical 
observations suggest a maximum of 2-8 days delay 1^ , 
and are compatible with no delay. 

Assuming the SN energy Esn ~ 2 x 10^^£'52.3 erg 
converts to kinetic energy with the observed velocity 
Vsnr ~ 3 X lO^wg.s cm/s of the SNR shell, the mass of 
the shell is Msnr « 2 — 3Mq. For our calculation, we 
take Msnr — 2.27712.2^0. The typical distance reached 
by the SNR shell is i?snr ~ lO^^ ^'wg.sid cm in td days. 
We assume an isotropic distribution of shell matter with 
d = Ai?siu./i?sni. = O.lS-i width for td > I day. For 
younger SNR shells, the matter may be taken to be ap- 
proximately isotropically distributed between the pro- 
genitor star and up to i?snr- 



One of the two dominant photon components trapped 
inside the SNR shell consists of photons from the SN 
shock with black body temperature of 0| 



^52. a*^-! ^ 



eV. 



(11) 



Another photon component may arise, for td ^ 1 day, 
because of a magnetohydrodynamic (MHD) wind which 
impacts the inner shell radius driving a forward shock 
through the shell. With an MHD wind luminosity of 
imhd 10^^im46 crg/s from a fast-rotating pulsar re- 
sulting from the SN, the thermalized photon energy in 
the SNR shell is 
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The optical depths for p7 interactions of these two photon 
components are 
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High energy protons escaping from GRB internal 
shocks (with an escape factor ?7p < 1) may interact with 
thermalized SN and MHD photons in the SNR shell. The 
corresponding neutrino threshold energies are 



Ef. 



El 



.th 



.th 



5.7 X lO'' 



1.7 X 10* 



pl/4 rl/3nl/4 
-^52.3"-! ^* 

3/4,1/2-p 
^9.5 '^rf I3OO 
.1/4.1/4 

-'^m46"-l 



(14) 



in the observer frame. Protons below A production 
threshold energies may undergo pp interactions in the 
SNR shell. The mean optical depth for such pp interac- 
tions is <Tpp>~ 60 Csh'm2.2l {vQ.5tdY average 60 
mb total pp cross section. Here C,sh ^ 1 is the fraction of 
cold protons inside the SNR shell. The maximum neu- 
trino energy from pion synchrotron energy losses in the 
SNR shell magnetic field is 
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We calculate the neutrino flux from pp interactions, 
below the proton threshold energy at A production, for 
times td^l day as 

,-2 VpCshfppfinLj^^iE^ 



E^^^ = 1.2 X 10 

"7 



-^27.2 



dE, 



M,{Ep) 
El 



GeV cm-2 s'^ 



M,{Ej,) 



Q. ( 1 ^ ^ 1 gp ^ 

7 ^ \^4 GcV /cm ^ ^ 4 GcV ^ 



V Gov / \^ Bp 



Ggv A 



(16) 



4 



where fpp^o « 1 is the neutrino conversion efficiency in 
pp interactions and E^"" « 20 x i;™'™ from Eq. (|Tl|l. 
The pj neutrino flux component is similar to the burst 
neutrino flux and we calculate it as 
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where /™q w 1 is the fraction of proton energy lost to 
pious. We assume that all protons (j^p = 1) below and 
10% protons {rjp ~ 0.1) above the proton threshold en- 
ergy 20 X [Eq. Q] escape the GRB internal shocks 
in Eqs. (|16|l and H17|) to calculate neutrino flux. 

For td ^ 1 day the calculation needs to include pion 
energy losses, due to inverse Compton scatterings in a 
high density photon field, before they decay to neutrinos. 



D. Precursor 

A precursor neutrino component (^^ 30 s duration) is 
generated while the GRB jet is still burrowing its way 
out of the stellar envelope. We have calculated this com- 
ponent for a similar progenitor model in detail elsewhere 
[l5j . There we considered a presupernova star which had 
either lost (i?* ^ 10^^ cm) or retained its H envelope 
(-R* ~ 10^^ cm). For the simultaneous SN-GRB sce- 
nario, these presupernova models are the same. In the 
"supranova" scenario, the status of the progenitor is un- 
certain, since numerical simulations which might provide 
guidance are not available so far. One possibility is that, 
except for the ejected envelope, the rest of the star be- 
comes the pulsar plus possibly a disk prior to collapse 
to a black hole leading to the GRB, so that a precursor 
may not be seen before the GRB event, since there is no 
stellar core to burrow through. A second possibility (de- 
pending on the delay) is that, at the time of the GRB, 
there is some fraction of the stellar gas which is not part 
of the ejected shell at radii intermediate between the shell 
and the central object, which could lead to a precursor 
resembling the above progenitor cases. The gas column 
density traversed by the jet is likely to be smaller than in 
the simultaneous SN case, reducing the pp contribution 
and producing a weaker precursor. 

We have plotted muon neutrino flux components at 
earth from GRB 030329 associated with SN 2003dh in 
Fig. ^ The "Precursor" components arrive ~ 30 s prior 
to the "Burst" component, and are expected in the case 
when the GRB and the SN events take place simultane- 
ously. The "Precursor I" signal assumes the presuper- 
nova stellar radius is « 10^^ cm. The "Precursor 11" 
signal is calculated assuming the star has an envelope, or 
the He/C/0 core has expanded, out to a radius ~ 10^^ 
cm before the onset of the GRB event. The "Afterglow 



(ISM)" and "Afterglow (wind)" components are calcu- 
lated when the GRB jet materials run into lower density 
ISM and a higher density wind respectively. While the 
burst and the afterglow components are based on generic 
GRB features, the "Supranova" components are highly 
model dependent. We have plotted the supranova com- 
ponents assuming that 10% of the shock accelerated pro- 
tons escape internal shocks to interact with the SNR shell 
which is up to 8 day old, according to current observa- 
tional limits. We have plotted the "Supranova" compo- 
nents for 0.1 d (short dashed curve), 1 d (long dashed 
curve) and 8 d (dotted curve). Note a reduced P7 flux 
contribution {E^, > 10^ GeV) for the 0.1 d case due to 
inverse Compton losses of pions. The duration of all com- 
ponents are At « 30Ai3o s. 
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FIG. 1: Muon neutrino flux components from GRB 030329 
and the associated SN 2003dh. Possible "Supranova" com- 
ponents are plotted for 0.1 d (short dashed curve), 1 d (long 
dashed curve) and 8 d (dotted curve) delays of the SN event 
prior to the GRB. The "Precursor" signals are calculated in 
the case when the GRB and the SN events take place si- 
multaneously. The "Precursor I" signal is for a progenitor 
star which has lost its H envelope. The "Precursor H" sig- 
nal assumes that the star has retained an H envelope, or the 
He/C/0 core has expanded, out to a radius ~ lO^'^ cm before 
the onset of GRB event. The internal shock "Burst" compo- 
nent [3 is shown by full dark hues. The "Afterglow (ISM)" 
and "Afterglow (wind)" components, also shown as full dark 
lines, are calculated for expansion into a lower density sur- 
rounding ISM, and a higher density wind, respectively. 



IV. DETECTION AT EARTH 

Neutrino detection at earth requires a large volume of 
material because of a small neutrino-nucleon (vN) cross 
section 21, 22, 23, 24J. Ice Cherenkov detectors such 
as AMANDA (251, IceCube 0, RICE [g] and ANITA 
0; water Cherenkov detector such as ANTARES 
and extended air shower detector such as AUGER 
are poised to detect cosmic neutrinos in different en- 
ergy ranges. Several authors have recently calculated 
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expected neutrino events, based on generic flux models, 
from individual GRBs |2y| . Here we calculate the number 
of neutrino events that one can expect from GRB 030329 
and the associated SN 2003dh, based on flux models de- 
scribed in Sec. IIIII in an under-ice Cherenkov detector 
located at the South Pole as an example. 

Neutrinos from a point source in the northern hemi- 
sphere travel through a substantial amount of earth ma- 
terial to reach a detector at the South Pole. Given a 
certain zenith angle 6, the column depth of material 
is Z{9) = p{r]9,l)dl, where p{r;9J) is the density 
of earth's interior, based on the Preliminary Reference 
Earth Model [23|, at a depth / and the corresponding 
earth's interior radius r. The upper limit of integration 
is the chord length L ~ 2Rq cos 9 where ~ 6371 km 
is earth's outer radius. The interaction length for neutri- 
nos inside earth is K^n{Ev) = Na<j'^%{E^), where Na is 
the Avogadro's number and the effective cross section is 
given by 
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for an incident neutrino flux <i>^. Here cr^^ and a 
are the total charge current (CC) and neutral current 
(NC) cross sections respectively and the third term cor- 
responds to the neutrino regeneration phenomenon due 
to NC interactions which degrade a neutrino to lower 
energies [23, H^- For simplicity we have ignored neu- 
trino regeneration in Eq. (|18|) . and the effective cross 
section is then the sum of CC and NC cross section 
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'^uN = '^ °N- The survival probability of 
a neutrino, from a point source, penetrating through the 
earth is 
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We use the total i/N cross sections (CC and NC) given 
in Ref. to calculate the survival probability. 

Neutrinos surviving after penetrating earth must inter- 
act with nucleons inside the effective detection volume to 
be detected as an event. Secondary e or /i produced by 
the CC interaction {ve,fi + N ^ e, p + X) then lead to 
either a high energy e-cascade or /Lt-track from which a 
Cherenkov si gnal can be picked up in the radio ^29] or 
optical band |25(. The probability for a I'e,^ to inter- 
act and produce a secondary e, p of energy Ee,^ above 
a minimum value i?™™ 
then 



in the vicinity of the detector is 
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Here Re is the electron range and is essentially the linear 
dimension of the detector. The muon range is given by 
the analytic formula 
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where a — 2 x 10^^ GeV cm^ g^^ and b 

113 ■ We have plotted neutrino interaction proba- 
bilities, in a km deep under-ice detector, after surviv- 
ing through earth [S{E^)PeAE,^)] from GRB 030329 
{9 = 68.4°) in Fig. [3 Because of their longer range, 
the probability for muon tracks (solid curves) are higher 
than for electron induced cascades (dashed curves). We 
have also plotted the probabilities for horizontal and up- 
ward neutrino events for comparison. 
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FIG. 2: Probability [S(Ev)Pe,fx{E,,)] of neutrino interaction 
in the detector volume, after penetrating through earth from 
GRB 030329, in a km deep under-ice detector at the South 
Pole. Also shown are the probabilities for upward and hor- 
izontal neutrinos for comparison. Probabilities for a muon 
event and an electron induced cascade event are shown by 
solid and dashed curves respectively. 

The number of neutrino events, given a flux — 
d^N^/dE^dt, in a detector of effective area A^s is then 



N, 



AeffAi / dE. 



d'N, 
dE^dt 



S{E,)P,AE.), (22) 



from Eqs. H19|) and H20|) . We have used a duration 
At — BOAiao s for all fluxes. We have tabulated (see 
Table the expected number of muon tracks and elec- 
tron cascades from GRB 030329 in an under-ice detector 
of Acs — km^ for different neutrino flux components. We 
used EfA = 100 GeV for these calculations. 



V. SUMMARY AND DISCUSSION 

We have calculated the > TeV neutrino signals from 
gamma-ray bursts resulting from the collapse of massive 
stars associated with supernovae, both in the case where 
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TABLE I: Neutrino events in a km scale under-ice detector 
at the South Pole from the GRB 030329 and the associated 
SN 2003dh. We have also calculated events assuming similar 
GRB-SN with declinations 90° (upward events denoted by |) 
and 0° (horizontal events denoted by —>) cases. We have not 
considered detector response for our calculation. 
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the supernova occurs simultaneously with, or hours to 
days in advance of, the 7-ray event. Longer SN-GRB 
time off-sets of vifeeks to months as considered in the 
"supranova" scenario which would have had inter- 
esting consequences for the interpretation of GRB after- 
glows, are ruled out by observations of GRB 030329-SN 



2003dh 0,0,0]. Shorter off-sets of days to hours, how- 
ever, which cannot be ruled out by optical observations, 
could be the signature of a delayed black hole formation 
from fall-back of gas onto an initial neutron star. Such 
two-step black hole formation could be expected in inter- 
mediate mass (3OM0 < < 401/©) stellar progenitor 
core collapse Q, after delays of fractions of hours to < 
day. Within this range, the time off-set is so far poorly 
constrained due to the approximate nature of the nu- 
merical simulations. Thus, any experimental constraints 
would be very useful, and neutrino signals such as dis- 
cussed here would be essentially the only way to obtain 
information about the dynamics of the core collapse. Fig- 
ure^gives the predicted supranova muon neutrino fluxes 
for a burst such as GRB 030329 assuming three different 
time delays between between the SN and the GRB. For 
a simultaneous GRB-SN event, these supranova fluxes 
would be absent. 

The detection of > TeV neutrino signals from tran- 
sient sources such as GRBs is basically background free 
0. The number of //-tracks and e-cascades from the 
GRB internal shock [Burst and Precursor II components] 
and from interaction of shock-accelerated GRB protons 
with the supernova shell (if the later precedes the GRB 
by hours to days) are certainly above atmospheric back- 
ground of ^ 0.3 events in the TeV-PeV range. In the 
PeV-EeV range, events from the Burst component 
should be detectable over zero atmospheric background 
events. We note, furthermore, that the number of 
events would be doubled in case there is no flavor oscil- 
lation. The detection of, or upper limits on, such events 
could provide important insights both on the collapse dy- 
namics and on the GRB central engine. 
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